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ABSTRACT 

A new procedure for smoothing a gamma-ray burst (GRB) lightcurve and calculating 
its variability is presented. Applying the procedure to a sample of 25 long GRBs, we 
have obtained a very tight correlation between the variability and the peak luminosity. 
The only significant outlier in the sample is GRB 030329. With this outlier excluded, 
the data scatter is reduced by a factor of ~ 3 compared to that of Guidorzi et al. 
(2005), measured by the deviation of fit. Possible causes for the outlier are discussed. 
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1 INTRODUCTION 

Despite many exciting progresses in observations, the nature 
of Gamma-Ray Bursts (GRB s) re mains to be a big puzzle 
[see Ivan Paradiis et all J2000) and IPirarJ ((2004) for recent 
reviews]. In such a situation, it is extremely important to 
identify some good correlations between the apparent (easy 
to measure or calculate) and intrinsic properties of GRBs. 
Several such correlations have indeed been found. For ex- 
ample, an anti-correlation between the the pe ak luminosity 
and sp ectral lag of GRBs has been found by iNorris et ID 
(2000), and a correlation between the peak luminosity and 
the variability of GRB lightcurves has been found by Feni- 
more & Ramirez-Ruiz (2000, hereafter FR00) and Reichart 
et al. (2001, hereafter R01). A correlation between the to- 
tal isotropic energy and the peak energy of the spectrum 
iAmati et al. 1 120021) . or the collimation-co rrected total en- 
ergy and the peak energy of the spectrum dGhirlanda et al] 
2004), has also been discovered. 

Recently, Guidorzi et al. (2005, hereafter G05) tested 
the correlation between the variability and peak luminosity 
of GRBs, using an expanded sample of 32 GRBs with mea- 
sured redshifts. The definitions of the variability and the 
peak luminosity are the same, but the size of the GRB sam- 
ple of G05 is about three times bigger than that of R01. The 
existence of a correlation was confirmed, but the scatter in 
the correlation is significantly larger than that found by R01 
(see, however, Reic hart fc Ny sewa nder 2005). Although the 
issue is in debate fGuidorzil 120051 : iReichart fc Nvsewanderl 
120051 : lReichardl2005t) . it is clear that with the definition of 
variability given by R01, the correlation between the vari- 
ability and the peak luminosity is not tight. 
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In this paper, we present a new definition for the vari- 
ability of GRB lightcurves. We then apply it to a sample of 
25 long duration GRBs with measured redshifts, whose data 
are publicly available. We show that, with the new definition 
of the variability, the correlation between the variability and 
the peak luminosity of GRBs is significantly improved: the 
data scatter is considerably reduced. 



2 NEW DEFINITION OF THE VARIABILITY 

To measure the variability of a lightcurve, first we must de- 
fine a reference lightcurve. The reference lightcurve should 
be sufficiently smoother than the original raw lightcurve. 
Since an ultimate model for GRBs does not exist yet, 
there is no first principle guiding us in choosing a reference 
lightcurve. What people usually do is to smooth the raw 
lightcurve with a linear "box car" filter (or moving window) , 
which smoothes the lightcurve with linear average (FR00; 
R01). 

Here we use a Savitzky-Golay filter iPress et al .120021) to 
smooth a lightcurve. The Savitzky-Golay filter is a more gen- 
eral and more powerful approach for smoothing noisy data 
than the linear box-car filter. The basic idea of Savitzky- 
Golay filtering is to approximate the underlying function 
(i.e., the reference lightcurve) within the moving window by 
a polynomial of higher order. An advantage of the Savitzky- 
Golay filter to the linear filter is that the former preserves 
high moments while the latter does not. 

A Savitzky-Golay filter is specified by three numbers: 
the order of the polynomial (m), the number of points used 
to the left of a data point (til), and the number of points 
used to the right of a data point (br) . To apply the Savitzky- 
Golay filter, the dat a must be binned w ith constant spacing. 
For more details see iPress et al.1 J2002I) . 
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We use a third order Savitzky-Golay filter. That is, we 
set m = 3. The width of the filter, i.e. np = til + np + 1, is 
determined with the approach of R01: np is set to be equal to 
the number of data points corresponding to a timescale Tf- 
the time spanned by the brightest 100/% of the total counts 
above the background (for details see R01). It turns out that 
/ = 0.5 most suits our purpose (rather than the / = 0.45 
used by R01 and G05; see Section 5). So, throughout the 
paper we use np determined by To. 5. 

We then define til = int[(np — l)/2], and np, = np — 
7jl — 1. If np is odd, we have np — til, i.e. the filter is 
symmetric about the point to smooth. If np is even, we 
have np = til + 1, the filter is asymmetric. In this case, 
we smooth the lightcurve twice: first we use the ul and np 
defined above, then we switch n^ and np. The results are 
then averaged. 

Suppose we have obtained a reference lightcurve by ap- 
plying the Savitzky-Golay filter to the raw lightcurve. Let 
us denote the count of the raw data in the i-th time bin by 
Ci, the count given by the reference lightcurve by Yi. The 
total squared deviation of the raw lightcurve from the refer- 
ence lightcurve is then J^-J^* 1 — ^) 2 j where ATbm is the 
total number of bins to be summed. We obtain the intrinsic 
squared deviation by subtracting the Poisson noise -/Vp i sson 

iVbin 

AC* 2 =W^2(d- Yif - TVpoisson , (1) 
i=l 

where W = np / (np — m — 1) is a statistical weight ac- 
counting for the fact that among the np data points only 
np — m — 1 are statistically independent. The inclusion of W 
allows us to apply the variability definition to any lightcurve 
with np > m + 1 . 

The summation in equation Q is from time t\ (corre- 
sponding to i = 1) to time £2 (corresponding to i = A^m), 
enclosing a major part of the lightcurve. Following FR00 
and R01, we define t\ to be the start of T90, ti to be the 
end of T90, where Tgo is the time during which the cumu- 
lative counts of the GRB increase fr om 5% to 95% above 
background iKouveliotou et al 

The Poisson noise is calculated by 

Whin 

iVpoisso* = \- Ci + (£ - ^CW.i] > ( 2 ) 

1=1 

where Cb g ,i is the background, the factor £ is the ratio of 
the background fluctuation to the Poisson noise of the back- 
ground (given by the reduced \ 2 °f the background fit). For 
GRBs detected by Swift we found that £ is often significantly 
larger than unity, indicating that the background fluctuation 
is quite non-Poissonian. 

The variability of the lightcurve is defined by the nor- 
malized squared deviation. We find that the following defini- 
tion leads to the tightest correlation between the variability 
and the peak luminosity 

AC 2 

V = (JVwn-l)C&« ' (3) 

where AC 2 / (N hin - 1) is the average of the squared devi- 
ation, C m ax is the net peak count (i.e., the background is 
subtracted). 

Our variability is defined in the observer's frame, so 



that the information of GRB redshift is not needed. This 
not only makes the computation of variability simple, but 
also eliminates an uncertainty arising from the assumption 
about the dependence of lightcurve variability on photon 
energy. The effect of GRB redshift was considered by R01 
and FR00 who defined their variabilities in the GRB frame, 
but the dependence of variability on redshift turned out to 
be very weak. 

As we applied our smoothing procedure to the GRBs 
in the sample described in the next section, we found that 
the tightest correlation between the variability and the peak 
luminosity is obtained if we iteratively apply the Savitzky- 
Golay filter Afiter times, where A/iter is the integer closest to 
Tgo/Tf (i.e., Alitor is roughly the number of moving windows 
contained in T90). 

In summary, our definition of variability differs from 
that of FR00 and R01 in the following aspects: (1) Our vari- 
ability is defined in the observer's frame, while the variabil- 
ities of FR00 and R01 are defined in the GRB's frame. (2) 
We normalize the average of the squared deviation by the 
squared peak count (the same as FR00), while R01 normal- 
ize the total squared deviation by the sum of squared counts. 
(3) FR00 and R01 use a linear box-car filter, while we use a 
nonlinear Savitzky-Golay filter. 



3 DESCRIPTION OF THE SAMPLE 

Our GRB sample contains 19 GRBs from the sample in 
G05, and 6 more GRBs detected by FREG ATE/ HETE-2 
and BAT/ Swift. So, the total number of GRBs in our sam- 
ple is 25. They are listed in Table with measured red- 
shift, calculated isotropic-equivalent peak luminosity, calcu- 
lated variability, and the number of iterations in applying 
the Savitzky-Golay filter. We have chosen to use the GRBs 
with data available publicly, which include GRBs detected 
by BATSE/CGRO, HETE-2, and BAT / 'Swift. 1 To obtain 
a reliable calculation of variability, we have only selected 
GRBs with more than 30% of total counts above the 3-cr of 
background. As a result, those GRBs with too low signal- 
to-noise ratios are not included in our sample. 

The 19 GRBs from G05 are: 970508, 971214, 980425, 
980703, 990123, 990506, 990510, 991216, 000131, 010921, 
020124, 020813, 030328, 030329, 041006, 050401, 050505, 
050525, and 050603. Their peak luminosities are taken from 
the same paper. The rest 13 GRBs in G05 are not included 
in our sample for various reasons: either their data are not 
publicly available (noticeably the 7 GRBs detected by Bep- 
poSAX but not by BATSE or HETE-2), or their data are 
incomplete or have too low signal-to-noise ratios. Although 
050315 and 050319 in G05 were detected by B AT / Swift, 
their data were not available to us since as this paper was 
written the archive of Swift only contained data taken after 
1 April 2005. 

The 6 newly added GRBs are 030115a, 030528, 050408 
(from HETE-2), 050802, 050803, and 050820a (from Swift). 
Their peak luminosities are calculated in exactly the same 

1 Although the lightcurve data of some GRBs detected by 
Konus/iyind are available from http://gcn.gsfc.nasa.gov/gcn/ 
konus_grbs.html, we have chosen not to use them since the data 
are highly incomplete in the part of pre- and post-burst. 
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Table 1. Variability versus peak luminosity for 25 GRBs with known redshift. 
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( a ) Measured redshift. 

< b > Mission: B (BATSE/ CGRO), H (HETE-2), BS (BeppoSAX), K (Konus/ WIND), BSw (BAT /Swift), U 
(Ulysses), S (SROSS-C), N (NEAR), R (RossiXTE), O (Mars Odyssey), RH (RHESSI): the data used are 
taken from the first mission mentioned. 

( c ) Number of iterations in applying the Savitzky-Golay filter (see Section 2). 

( d ) Calculated variability (equation [3} and error. 

( e ) Isotropic-equivalent peak luminosity in 10 50 erg s — 1 in the rest-frame 100-1000 keV band, for peak fluxes 
measured on a 1-s time-scale (Hq = 65 km s" 1 Mpc -1 , Q m = 0.3, and Q/± = 0.7). 

W Reference s for t he r edshift measurem e nts: ( 1) iMetzger et alJ <1997tl. ( 21 IKulkarni et alJ Jl998f), (3) 
[ Tinnev et all 1199^1. (4) IDiorgovski et alJ tl998l). (51 IKulkarni et alJ <199^, (6) iBloom et alJ i2003h. (71 
beuermann et alJ jl999ll (81 IVreeswiik et alJ il999ft. ( 91 lAndersen et alJ I200CTI . ClOl iDiorgovski et alJ 

l200lll . (llllHiorth et alJ i2003ari. (12l|Price et alJ i2002t>. (131 http:// s pace.m it.ed u/HETE/Bursts/ Data/ 

andlSmith et alJ ^2005^ . (1 41 iMartini et alJ i2003l) , (1 51 iGreiner et all i2003l1 . (1 6llRau et alJ <2q05l). (17 ) 
iFugazza et alJ i2004ll (1 8) iFvnbo et alJ i2005a>). (191 iBerger et alJ <200^ an d iProchaska et alJ j20053l . 

(20l| Berger et al.||2005H1 . (21 1lFolev et al.U2005ll, (22 l lBerger fc BeckerH200,^ (23l lFvnbo et ali l2005bll . 

(24l lBloom et alJ 1200^ (251 IProchaska et alJ <2005bl) and lLedoux et alJ 120031 . 



way as in R01 and G05. For GRB s detected by HETE-2 , 
we used the spectral fit provided bv lSakamoto et, all (1200.41 
and the HETE website. For GRBs detected by Swift, peak 
spectra are extracted from their event files, then fitted with 
a power law which is sufficient for all cases (in the energy 
range 15-350 keV). 

Althoug h iGolenetskii et all feOOEl) . iPal'shin etafl 
l|2005l) . and ICumming^T^n ll^OOa) have reported the 
detection of a second, larger episode of emission from GRB 
050820a (by both Konus/ Wind and BAT /Swift), only the 
data for the first episode from Swift are available to us. 
Thus, for GRB 050820a, the variability and peak luminosity 
listed in Table 1 refer to its first episode (of duration ~ 30 
sec). 

The time bin of each GRB lightcurve is as follow. For 
GRBs detected by BATSE, except 000131, the time bin is 
64-ms. GRB 000131 does not have a 64-ms lightcurve, so 



we use its 1.024-sec lightcurve. Although a bin of 1.024-sec 
is somewhat large, the duration of 000131 is about 100-sec 
which is much larger than its time bin. For GRBs detected 
by HETE-2, the time bin is 164-ms. While for GRBs de- 
tected by BAT /Swift, we extract 64-ms lightcurves from 
their event files available in the public data archive. 

As we analyze the GRB lightcurves we use the total 
counts in the energy range specified as follow: 25-300 keV 
for GRBs detected by BATSE, 30-400 keV for those de- 
tected by HETE-2, and 25-350 keV for those detected by 
BAT/ Swift. We have tried to make the range of energy for 
GRBs detected by different instruments have largest over- 
lap. 

Finally, we remark that for the GRBs detected by 
BATSE, the lightcurve data in the pre- and post-burst sec- 
tions were divided equally into 64-ms time bins from their 
original 1.024-sec time bins. When we calculate the Poisson 
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Figure 1. Peak luminosity versus variability of 25 GRBs. The solid straight line is the least-x 2 linear fit to the data, with GRBs 980425, 
030329, and 030528 excluded (see the text): log L = 3.25 log V + 59.42 (L in units of erg/sec). The x 2 /dof = 38.58/20=1.93. The two 
dashed lines mark the 1-<t deviation of the fit. (Filled circles: GRBs detected by BATSE; Squares: GRBs detected by HETE-2; Triangles: 
GRBs detected by BAT /Swift.) 



noise in those sections, we multiply the result by a factor of 
0.0625 (= 0.064/1.024) to take into account the reduction 
in the Poisson noise arising from the change in time bin. 



4 RESULTS 

We have applied the smoothing procedure described in Sec- 
tion 2 to our GRB sample. The obtained variabilities and 
errors, as well as the number of iterations in applying the 
Savitzky-Golay filter, are listed in Tabled The errors of vari- 
ability are principally caused by photon noises. However, the 
uncertainties arising from changing .Niter (see Section 2) by 
±1 are also taken into account: we find the changes in the 
variability for Niter ± 1 an d add in quadrature the maximum 
to the statistical error. 

The peak luminosity versus variability is shown in 
Fig.Q Clearly a very tight correlation between the two quan- 
tities exists, with only one prominent outlier: GRB 030329. 
Given the large error in its variability, it is unclear if GRB 



980425 — which is famous for its smallest redshift, lowest 
peak lumino sity and least total en ergy, and association with 
SN 1998bw jGalama et al]ll99gft — is off trend as in other 
kind of correlations (e.g., the total energy-peak energy cor- 
relation). In fact, the error in variability is larger than or 
about the same value as the variability itself for 980425 and 
030528 (see Table 1). 

We have made a least-x 2 linear fit to log L— log V, where 
L is the peak luminosity in erg/sec. To take into account 
both th e errors in L a nd V , we have made use of the program 
fitexy in lPress et alJ i2002l) . The asymmetric errors in logL 
and log V are averaged. That is, we take 



[log (q + a q ) - log (q - a q )] 



(4) 



where q = L or V. Since the errors in log V are infinite 
for 980425 and 030528, these two GRBs do not contribute 
to the total x 2 - GRB 030329 was not considered due to its 
large offset. With 980425, 030329, and 030528 excluded, the 
total number of GRBs in the final sample to fit is 22. We 
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Figure 2. The confidence intervals in the a — b space. The cross 
sign denotes the values of a and b at the minimum of x 2 (a = 3.25, 
b = 59.42, x 2 = 38.58). The contours correspond to the increase in 
X 2 , Ax 2 = 2.30, 4.61, and 6.17 respectively. The region enclosed 
by each contour contains 68.3%, 90%, and 95.4% of normally 
distributed data. 

have obtained the following result (the solid straight line in 
Fig. ID 

logL = alogV + &, (5) 
where 

a = 3.25 ±0.26, b = 59.42 ± 0.53 . (6) 

The reduced chi-square \r = xVdof = 1.93, where the de- 
gree of freedom dof = 22 — 2 = 20. The smallness of Xr 
indicates a very tight correlation between the peak luminos- 
ity and variability, for details see Section 5. 

The two dashed lines in Fig. mark the 1-a width of 
the fit, measuring the scatter of data and denned by the 
deviation of fit (see Section 5). However, this should not be 
mixed with the error in the predicted log L, which we discuss 
below. 

In Fig. [5] we plot the confidence intervals in the a — b 
plane. The best fit values of a and 6, which correspond to 
the minimum of x 2 > are indicated by the cross sign in the 
figure. The three contours enclose the regions containing, 
respectively, 68.3%, 90%, and 95.4% of normally distributed 
data. The highly elongated shape of the contours indicates 
that a and b are highly correlated. 

For any measured value of log V, equation predicts 
a value of log L. The error in the predicted log L comes from 
two sources: the error in log V, and the uncertainties in a 
and b 

CT log L = ^ 2 <V+' T a6 • ( 7 ) 

Generally, the errors are asymmetric in the upward and 
downward directions. Let us denote by a a b,i the upward 
error due to uncertainties in a and b, by a a b,i the downward 
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log Luminosity (erg/sec) 

Figure 3. The error in the log L predicted by equation JSJi due to 
the uncertainties in a and b. Filled circles are numerical solutions 
for the upward error (a a b i). Open circles are numerical solutions 
for the downward error (cr a b,2)- Solid lines are analytical approxi- 
mations to the numerical solutions given by equations IAH - IXoV 

error. That is, the predicted log luminosity is logL^ a ''*, 
if <J\ og y = 0. In Fig. [3] we show the numerically calculated 
values of <T a i>.i (filled circles) and a a b t 2 (open circles) versus 
the predicted log luminosity. From Fig. [5] we find that 5b ~ 
2.05 5a when 5a and 5b are not close to zero, which leads to 
&ab ~ (log V + 2.05)5a. Thus, a a b is minimized at log V m 
— 2.05, i.e. at logi ~ 52.8, consistent with the numerical 
results. 

For the convenience of application of the log L — log V 
relation JSJ to the prediction of GRB's peak luminosity, we 
have attempted to approximate the numerical results for a a b 
by several analytic formulae. The results are summarized in 
Appendix A, and represented by the solid lines in Fig. 

Now we apply the above results to GRB 030329 to check 
how far it deviates from the log L — log V relation. From Ta- 
ble 1, for GRB 030329 we have V = 0.0070 and a v = 0.0010. 
Submitting \ogV — —2.15 into equation we obtain a 
predicted log luminosity for GRB 030329: logL = 52.4. 
Since GRB 030329 lies below the straight line in Fig. □ the 
relevant error in the predicted log luminosity is determined 
by Oab,2 and a\ og v = logV — log (V — av) « 0.067. By 
equation 1A4II we have o" a 6,2 ~ 0.08. Then, by equation J7J, 
the total error in the predicted luminosity is ai og i, ~ 0.23. 
From Table 1, the measured luminosity of GRB 030329 is 
logL « 50.8, with an upward measured error <7i og i, ~ 0.08 
that is smaller than the error in the predicted log luminosity. 
So, GRB 030329 deviates from the correlation (the straight 
line in Fig.[TJ by (52.4 - 50.8)/0.23 « 7a. 



5 COMPARISON TO GUIDORZI ET AL.'S 
RESULTS 

To make a quantitative comparison of our results to that of 
G05, we need a parameter that measures the scatter of data 
points around the best fit model. The parameter should not 
be sensitive to the variat ion in data error s caused by different 
definitions of variability. iReichard (1200 ll) has provided such a 
parameter based on the maximum-likelihood method, which 
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Figure 4. Definition of the variance of fit. The thick solid curve 
is the fit y = f(x). The j/-component of the variance of fit, s y , 
is defined by equation 1111 , The two dashed curves, called "the 
1-<t width" of y = f(x), are defined by y = f(x) ± s y . The x- 
component of the variance of fit, s^., is related to Sy by equa- 
tion 1141 . The scatter of data points around y = f(x) is measured 
by s = s x s y / + s y , the length of the straight line OC that 
is perpendicular to the straight line AB. [The straight line AB 
locally approximates the curve y = f(x) when s y is sufficiently 
small or the curvature of the curve is negligible.] 



is called the "sample variance" , or "slop" . However, here we 
choose to u se the "variance of fit" d efined with the approach 
of least-x 2 feevington fc Robinsonll992l) . which is much less 
sophisticated but enough for our current purpose. 

Suppose a data set {xi,yi}, with errors {<J x ,i, fy,i} , is 
fitted by a model y — f(x). The fit minimizes the x 2 denned 
by 



AT 
i — 1 1 



where N is the total number of data points, and 

df^ 2 



2 2 
G% = a v.i + 



2 

x,i 



The reduced x 2 > i- e -? the x 2 per degree of freedom, is 
X l =X 2 /dof = X 2 /(N-P) , 



(8) 



(9) 



(10) 



where p is the number of parameters in the function f(x) 

The vari ance of the fit, measured along the y-direction, 
is defined by feeving"ton fc Robinsonlll992 



y ~ N - p ^ af 



-slVi- f {Xi)} , 



(11) 



JV 



(12) 



The variance of fit measures the spread of data around the 
model. We denote the variance of fit by s 2 (then s is the 
deviation of fit), to distinguish it from the data variance a 2 
which represents the error in the measurement of each data 
point. 

The appearance of (a 2 ) in the numerator and of in the 
denominator in equation (I11H results that the variance of fit 
is not sensitive to the variation in data variance. 

By the definition of Xr i equation 111 It can be recast into 
a simpler form 



2 / Z\ z 



(13) 



We define the variance of the fit measured along the 
x-direction by 

t -*(£)". cm 

where df/dx is the slope of the curve y = f(x). Note, by 
definitions, s y is always a constant, but s x is a function of 
x unless f(x) is a linear function of x. 

We measure the scatter of data points around the curve 
y = f(x) by the variance of fit along a direction locally 
perpendicular to the curve, which is (see Fig. 



2 2 

S x Sy 



1 + (df/dx) 



(15) 



For the problem in this paper, we have x — log V , y — 
log L, and f(x) — ax + b. Thus, df/dx = a is a constant, and 
s x = Sy/a, a i — a yi + a a x ^. The deviation of fit is then 

(16) 



where (of) is the weighted average of the individual data 



With the above tools, we are ready to compare our re- 
sults to that of G05 quantitatively. To make the comparison 
well denned, we analyze the overlapping sample which con- 
tains 19 GRBs (see Section 3). In Fig. we plot the vari- 
ability with our definition versus that with ROl's definition 
(the latter was copied from G05). The two variabilities are 
highly correlated, but the scatter is large. This large scatter 
is necessary for obtaining a better correlation between the 
peak luminosity and the variability of GRBs. 

The largest scatter in Fig. 03 occurs at GRB 980425, 
for which our background fit leads to T90 = 60.29 sec, in 
contrast to the 34.88 sec reported in the BATSE web page. 
If we adopted T90 = 34.88 and £ = 1 (£ is the ratio of the 
background fluctuation to the Poisson noise, see Section 2), 
we would obtain V = 0.00053 ± 0.00057 for 980425. This 
value of V is only slightly larger than that listed in Table 1, 
so 980425 remains to have the largest scatter in Fig. 03 With 
our definition of variability, GRB 980425 is much less offset 
from the variability-luminosity correlation. 

Figure 03 also shows that the values of the variability 
with our definition are systematically smaller than that with 
ROl's definition, as indicated by the straight line denned by 
"our variability = ROl's variability divided by 15". This is 
mainly caused by the different normalization in the two def- 
initions. In addition, our variability generally has a larger 
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Figure 5. The variability with our definition versus the variabil- 
ity with Reichart et al.'s definition for 19 GRBs. The straight 
line is defined by the equation "our variability = Reichart et al.'s 
variability divided by 15". 



error bar than that of R01, caused by the fact that we nor- 
malize our variability by peak count, which suffers larger 
photon noise than total counts; and we have considered the 
error arising from changing TVitcr by ±1. 

Now, we make a linear fit to log L — log V for the over- 
lapping sample, with our definition of variability. For the 
reasons explained in Section 4, we exclude GRBs 980425 
and 030329 during the fit. Thus, the total number of GRBs 
to fit is TV = 17. We obtain logL = 3.10 log V + 59.13, and 
X 2 /dof = 2.43 where dof = 15. The weighted average of 
the individual data variance is (of) = 0.02715. Thus, for 
the fit with our definition of variability, we have s y — 0.257, 
s x = 0.0829, and s = 0.079. 

Next, we make a linear fit to log L — log V for the over- 
lapping sample, with ROl's definition of variability. Again, 
GRBs 980425 and 030329 are excluded, to make the GRB 
members in the sample remain the same. We obtain logL = 
1.77 log V + 54.06, and xVdof = 19.89. The weighted av- 
erage of the individual data variance is (of) = 0.009240. 
Thus, for the fit with ROl's definition of variability, we have 
s y = 0.429, s x = 0.2417, and s = 0.211. 

The scatter of data points around the fitted log log V 
is measured by s. Thus, according to the above numbers, for 
the overlapping sample the scatter in our log L — log V rela- 
tion is smaller than that of G05 by a factor of 0.211/0.079 = 
2.7, although the weighted average of the individual data 
error (= (of ) 1 ^ 2 ) with our definition of variability is larger 
than that with ROl's definition of variability by a factor of 
(0.02715/0.009240) 1/2 = 1.7. 

The above results, for the overlapping sample excluding 
980425 and 030329, are summarized in Fig. |SJ where filled 
circles and thick lines represent the results with our data, 
open circles and thin lines represent the results with G05's 
data. The luminosities in the two set of data are the same, 
so each GRB is represented by a horizontal line in Fig. |S| 
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Figure 6. Comparison of our results (filled circles and thick lines) 
to that of G05 (open circles and thin lines). Straight lines repre- 
sent the least-x 2 fit and the 1-cr width. (Our variability has been 
multiplied by a factor 15.) 
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Figure 7. Variation of the deviation of the linear fit to log L — 
log V, with respect to the choice of N- ltcI . The parameter a, which 
the To. 5 in the definition of iVi ter is multiplied by, is set to be the 
same for all GRBs in the sample. The variabilities presented in 
this paper correspond to the case a = 1. 



However, the variabilities in the two set of data differ, so we 
see the offset along the horizontal direction of each pair of 
filled circle-open circle. The solid straight lines are least-x 2 
fits to each data set, and the dashed straight lines mark the 
"l-<7 width" of the fits. We have set the unitary length of the 
horizontal and vertical axes to be the same, so that for each 
data set s is just the half-width of the region bounded by 
the two dashed lines that you see. The figure clearly shows 
that the scatter of our data set is much narrower than that 
of G05's data set. 

With the above defined scatter parameter s, we can 
also check how sensitive our results are to the parameters in 
our definition of variability. As mentioned in Section 2, we 
have chosen T/ = To. 5, and Niter to be the integer closest to 
the ratio Tgo/Tf. With this choice, for the GRB sample de- 
scribed in Section 3 (excluding 980425, 030329, and 030528), 
we obtained x 2 /dof = 1.93 and s = 0.078 (Fig. Q where the 
two dashed lines are obtained by shifting the solid line up- 
ward/downward by s y = 0.265). If we followed R01 and G04 
to choose Tf = To. 45, we would obtain x 2 /dof = 2.53 and 
s = 0.087 in the linear fit to log L — log V . A slightly larger 
value of s indicates that the correlation is slightly looser. 

To test how sensitive the results are to the choice of 
•Niter , let us redefine .Niter by the integer closest to T90 / alb. 5 , 
where a is a parameter that is set to be the same for all 
GRBs in the sample. When a = 1, we return to the definition 
that we have adopted. We now let a vary from 0.2 to 2, and 
calculate the corresponding s in the least-x 2 fit to logT — 
log V. The calculated s versus a are presented in Fig. |7| 
which shows that the results are stable to the choice of -Niter- 
That is, the variance of fit does not change drastically with 
a small change in a. The figure also shows that a = 1 is 
around the minimum of s(a). 



6 CONCLUSIONS AND DISCUSSION 

We have presented a new definition for the variability of 
GRB lightcurves. The new variability is defined in the ob- 
server's frame, and the lightcurve is smoothed with the 
Savitzky-Golay filter. The former makes it easy to apply 
the procedure to all GRBs with well measured lightcurves 
(i.e., with a significant signal-to-noise ratio, and a reasonable 
size of time bin — neither too large to erase the variability on 
short timescales, nor too small to introduce too large Poisson 
noises), without referring to the redshift information. 

We have applied our smoothing procedure to a sam- 
ple of 25 long duration GRBs with measured redshifts and 
publicly available data. A very tight correlation between the 
variability and the peak luminosity is found, with only one 
prominent outlier: GRB 030329. Excluding 030329, 980425, 
and 030528 (the latter two are due to their too large errors 
in variability), a linear fit to the log of variability and the 
log of peak luminosity is obtained (equations |S] and with 
X 2 /dof = 1.93. The smallness of the reduced \ 2 indicates 
that the data scatter is very small (Fig. [T]and Fig.®. 

We note that, although the existence of a sec- 
ond, larger e pisode of emission from GRB 050820 a has 
been claimed jGolenetskii et alJl2005t iPal'shin et"afl 120051: 
ICummings et alH2005li ~and we have only made use of the 
data for its first episode, GRB 050820a fits the relation per- 
fectly well. 

Our results are a significant improvement to that of 
G05 (as well as that of FR00 and R01), who used a different 
definition of variability and obtained xVdof = 1167/30 = 
38.9. An analysis on the overlapping sample of 17 GRBs 
(after GRBs 980425 and 030329 being excluded) shows that 
the scatter in our data is smaller than that in G05's data by 
a factor of 2.7, measured by the deviation of fit (Section 5). 

The improvement to the correlation is caused by using 
a new smoothing procedure and a different normalization in 
the definition of our variability, not by neglecting the effect 
of GRB redshift. R01 have shown that the dependence of 
variability on the redshift is extremely weak. The fact that 
a very tight correlation between the peak luminosity and the 
variability is obtained without including the corrections to 
the variability from the GRB redshift also leads us to believe 
that such corrections are not necessary. 

A remarkable feature of the correlation that we have 
found is that it does not rely on the corrections to the lumi- 
nosity from the collimation of GRB jets. In other words, 
the peak luminosity used above is simply the isotropic- 
equivalent peak luminosity. If this correlation is confirmed 
by future bursts, it will provide a convenient calibration to 
GRB luminosity and distance. 

GRB 030329, which is considered to b e a firm case 
for the connection of GR Bs with supernovae iStanek et alJ 
2003; Hior th et alJl2003bf) . is offset from the correlation by 

7 a (Fig. 1). Although the cause is not clear, we have no- 
ticed that 030329 and 050525 differ from other GRBs in the 
sample in the following way: their lightcurves consist of two 
distinctly separated pulses with each containing at least 30% 
of total counts. Motivated by this observation and the fact 
that GRB 050820a fits the relation well although only the 
first episode of its emission has been used, we have recal- 
culated the variabilities of 030329 and 050525 by dividing 
the smoothing window by a factor of 2 (i.e., dividing the 
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np obtained from To. 5 by 2) but keeping iViter unchanged. 2 
We obtained V = 0.0022 ± 0.0003 for 030329, which well 
fits the solid straight line in Fig. 1. For 050525, we obtained 
V = 0.0032 ± 0.0007, which fits the straight line equally 
well as the value listed in Table f . Whether this treatment 
is correct must be tested when more GRBs with lightcurves 
similar to that of 030329 and 050525 are available. 
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APPENDIX A: ANALYTIC APPROXIMATION 
TO THE ERROR IN THE LOG LUMINOSITY 

The numerical results for the error in the log L predicted by 
equation JSJ due to the uncertainties in a and b, a a t, can 
be approximated by the following formulae (a a b,i for the 
upward error, a a t,2 for the downward error): 
(1) When log L < 52 

<r aM « 7.2017 - 0.16777a; + 0.0006006a; 2 , (Al) 
o- ab , 2 « 8.0929 - 0.18000x + 0.0005134a; 2 , (A2) 

where x = log L, L is in units of erg/sec. (2) When 52 ^ 
log L < 54 

cr aM » 2274.034248 - 124.984716a; + 2.286986a; 2 

-0.013930421a; 3 , (A3) 

a abt2 » 3028.186992 - 167.524125a; + 3.087093a; 2 

-0.018948436a; 3 . (A4) 

(3) When log L > 54 

T a b,l ~ -6.8502 + 0.12989a; , (A5) 
Tab,2 « -5.7594 + 0.10933a; . (A6) 

The fractional errors in a a b given by these formulae are 
always < 5%. 

This paper has been typeset from a TjtX/ DTffX file prepared 
by the author. 



2 Roughly speaking, this is equivalent to smooth each GRB pulse 
separately. 



